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Abstract 

For the Lorentz gauge the influence of various Gribov gauge copies on the 
fermion propagator is investigated in quenched compact lattice QED. Within 
the Coulomb phase besides double Dirac sheets the zero-momentum modes of 
the gauge fields are shown to cause the propagator to deviate strongly from 
the perturbatively expected behaviour. The standard way to extract the fermion 
mass fails. The recently proposed zero-momentum Lorentz gauge is demonstrated 
to cure the problem. 

1 Introduction 



Lattice gauge theories allow to compute most of the relevant observables without 
any gauge fixing. Nevertheless, computations of gauge dependent objects, e.g. 
gauge or fermion correlators, can give us more detailed information about non- 
perturbative properties of quantum fields and allow a direct comparison with the 
perturbation theory in the continuum. 

However, it is well known that gauge fixing, in particular the Lorentz (or 
Landau) gauge, leads to the occurence of gauge or Gribov copies [0]]. For QED 
this happens even in the continuum, as long as the theory is defined with toroidal 
boundary conditions 0. 

In this paper we are going to consider compact U(l) lattice gauge theory with 
Wilson fermions || in the quenched approximation (qQED). In the weak coupling 
region, i.e. within the Coulomb phase, it describes massless photons (weakly) 
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interacting with fermions. Full lattice QED is expected to reproduce continuum 
perturbative QED, which has been experimentally proven with greatest precision. 

The standard iterative way to fix the Lorentz gauge for compact U(l) lattice 
gauge theory has been shown to lead to serious Gribov copy effects [§], |5|, [| [7|]. 
As a consequence the transverse non-zero momentum photon correlator does not 
reproduce the perturbatively expected zero-mass behaviour. For the fermion cor- 
relator a strong dependence on the achieved gauge copies has been reported, too 
[Q. The standard fermion mass determination becomes badly defined. Careful 
numerical || |7|, || f| and analytical [|10[ studies have shown that the main 
gauge field excitations responsible for the occurence of disturbing gauge copies 
are double Dirac sheets (DDS) and zero-momentum modes (ZMM). 

For the aim of achieving a better agreement between numerical lattice results 
and lattice perturbation theory, the Lorentz gauge fixing procedure was coupled 
with additional gauge fixing conditions. The authors of || proposed to employ 
a unique initial gauge realized by the "maximal tree" axial gauge condition. In 
|j] a non-periodic gauge determined from spatial Polyakov loop averages was 
proposed to suppress unwanted DDS. It turned out that a gauge which allows 
to suppress DDS is sufficient to produce a correct non-zero momentum photon 
correlator perfectly compatible with a vanishing photon mass. However, in this 
way the Gribov problem is not yet solved, because the removal of DDS does not 
automatically lead to the global extremum of the Lorentz gauge functional. As 
we have seen recently |§, the above mentioned axial gauge even fails to remove 
DDS. In the same paper we have shown that in order to solve the Gribov problem 
the ZMM of the gauge fields have to be necessarily suppressed, too. An alter- 
nating combination of the Lorentz gauge fixing steps with non-periodic gauge 
transformations suppressing ZMM provides a practical solution of the problem. 

In the present paper, we are mainly studying the influence of the ZMM on 
the fermion correlator when applying the Lorentz gauge. The aim is to check 
whether the results of |6|, 0, §, 0, [31], [□]] apply to the fermion correlator, too. 
DDS will be removed from the beginning in order to concentrate on the effect of 
the ZMM. We shall demonstrate that only the correct account of the ZMM will 
allow us to determine the fermion mass properly. 

The outline of our paper is as follows. In Section 2 we introduce our notations 
in particular for the fermion propagator. For the latter we provide the analytical 
representation in a fixed constant gauge field background. Section 3 presents the 
details of the gauge fixing procedures employed. In Section 4 Monte Carlo results 
for the fermion correlator obtained with different Lorentz gauge fixing procedures 
are compared with the zero-momentum mode approximation of the correlator, for 
which only the constant background gauge field modes are taken into account. 
The conclusions are drawn in Section 5. 
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2 The Fermion Correlator 



We consider 4d compact U(l) gauge theory on a finite lattice with size V = 
iV s 3 x N t . The standard Wilson action consists of the pure gauge Sq and the 
fermion contribution Sp as follows 

S G = f3 E (1-cos^). (2.1) 

X,[l<U 

Ox,w = &x,ii + Ox+n, v - Qx+o,» - 9x,v is the plaquette angle with 6 X ^ E (-vr, vr] 
denoting the link gauge field variable, f3 = I /el is the inverse coupling. The 
lattice spacing is put a = 1. 
The fermion part is given by 

S F = J2^ M xyWy (2.2) 

= (4 + m ) E^A - \ E{^e if Ml - +^+A e " i9 " M ( 1 + 

X X,[l 

ip x denoting the fermion field, 7^ the Hermitian Dirac matrices. 

For the gauge and the fermion field we apply periodic boundary conditions 
(b.c.) except for the fermion field in the imaginary time direction x 4 , where we 
prefer to use anti-periodic ones. 

We are going to study the fermion correlator for given gauge fields 9 Xtfi 

r(r; 6) = 1 E E M *~, 1/4 = *4 + r. (2.3) 

x,X4 y 

For simplicity, we restrict ourselves to the scalar and vector parts of the fermion 
correlator, respectively 



T s (r;9) = ^ReTr (r(r;0)); 

r y (r;^) = ^ReTr ( 74 r(r^)), (2.4) 

where the trace is taken with respect to the spinor indices. For anti-periodic b.c. 
in x 4 the vector (scalar) part becomes an even (odd) function in r around 
r = Nt/2 , for periodic b.c. vice versa. 

In qQED the above correlator has to be averaged with respect to the gauge 
field 9 X ^ with the weight exp(— Sq) • Lateron, we shall compare the quan- 
tum average ( T )g with the zero-momentum mode approximation where only 
background gauge fields being constant in space-time are taken into account. 
Therefore, we construct analytically the correlator for a (uniform) gauge config- 
uration given by 

Ox,n = 4>^ -7T < M < 7T, fj, = 1, • • • , 4. (2.5) 
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One obtains the following finite size results for the scalar and vector parts, 
respectively 



r.cfr;. 



x 



Jt,0 



l + £ 2 -2£(1 + M) 



x 



2(1 + M) l-£ 2 
[£ T - £ 2N ^] cos(0 4 r) - c[£ Nt+T - £ Nt - T ] cos[<p 4 (N t - r)] 
1 + £ 2N > - 2c£ N * cos((f> 4 N t ) 



(2.6) 



rVfr; 



x 



1 - 5 Tfi 
2(1 + M) X 

[£ T + g 2JVt ~ T ] cos(0 4 r) - c[£ Nt+T + S Nt ~ T ] cos[0 4 (iV f - r)] 
1 + £ 2N * - 2c£ N * cos((j) 4 N t ) 



(2.7) 



where c = — 1 (c = +l) holds for antiperiodic (periodic) boundary conditions 
in x 4 and 



£ = l + 



M 2 + IC 2 VM 2 + JC 2 ^/(M + 2) 2 + lC 
2(1 + M) 



2(1 + M) 



M. = m + (1 — cos0/) 



i=i 



\ i=i 



sin 



m > 0. 



If we put all 



0, the Eqs. ( |2.6| , |277|) reproduce the results for the free 



fermion correlator [12 



In practical computations the fermion field ip x is rescaled with a factor 
1/(4 + mo) 1 ' 2 and the bare mass value mo is replaced by the hopping-parameter 
k by 

k=— — -. (2.8) 

2(4 + m ) 1 ; 



3 Lorentz Gauge Fixing 



A lattice discretization of the Lorentz (or Landau) gauge fixing condition can be 
written as follows 



(si n ~ sin 6» x ._ /t>M ) = 0. 



(3.1) 



In practice, instead of solving this local condition one maximizes iteratively the 
gauge functional 

(3.2) 



F[9] = —Y,^s9 x ^ = M a x. 
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with respect to (periodic) gauge transformations 



Ox,n — * ®x,u = a x + Ox,n ~ a x+fL mod 2vr, a x € (-7T, vr]. (3.3) 



In our simulations the maximization of the gauge functional ( |3.2| ) has been 
continued until both the mean and the local maximal absolute values of the 
l.h.s. in eq.( |3.1| ) became less than some small given numbers (in our case, 1CT 6 
and 10~ 5 , respectively). In order to accelerate the maximization one can apply 
overrelaxation optimized with respect to some parameter u [ D| . 



We call the algorithm standard Lorentz gauge fixing, if it consists only of local 
maximization and overrelaxation steps. It is well-known that this procedure 
normally gets stuck into local maxima of the gauge functional. The solutions 
corresponding to different local maxima are called Gribov or gauge copies. 

It is a common believe (see also [H) that the Gribov problem has to be solved 
by searching for the global maximum of the gauge functional (|3.2|) providing the 
best gauge copy (or copies, in case of degeneracy). In || we have shown that in 
order to reach the global maximum we have necessarily to remove both the DDS 
and the ZMM from the gauge fields. 

DDS are identified as follows. The plaquette angle fl2.1|) can be decomposed 
[PH 6 X)I1V = ~Q x ^v + 27m XjAtl , , where the gauge invariant 6 X ^ U e (— 7T, 7r] can be 
interpreted as physical (electro-) magnetic flux and the discrete gauge dependent 
contribution 2irn x ^ u , n x ^ v = 0, ±1, ±2 represents a Dirac string passing through 
the given plaquette if n x ^ u ^ (the Dirac plaquette). A set of Dirac plaquettes 
providing a world sheet of a Dirac string on the dual lattice is called Dirac sheet. 
Double Dirac sheets (DDS) consist of two sheets with opposite flux orientation 
which cover the whole lattice and close by the periodic boundary conditions. 
Thus, they can easily be identified by counting the total number of Dirac pla- 
quettes Npp for every plane (/i; v). The necessary condition for the appearance 
of a DDS is that at least for one of the six planes (/i; v) holds 

<;> > 2^. (3.4) 

DDS can be removed by periodic gauge transformations. But - as it was demon- 
strated in H - the standard Lorentz gauge fixing procedure usually does not 
succeed in doing this. DDS occur quite independently of the lattice size and 
the chosen f3. As a consequence the numerical result for the non-zero momen- 
tum transverse photon correlator significantly differs from the expected zero-mass 
perturbative propagator [p], p], |10 . 



The ZMM of the gauge field 



* X 
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do not contribute to the gauge field action ( [2.1| ) either. For gauge configurations 
representing a small fluctuation around constant modes 

9 x , fJ , = 4>fi + 80 x ,n, = 0, \89x,n\ < 1, 

X 

it is easy to see, that the maximum of the functional Q3.2|) requires = . The 
latter condition can be achieved by non-periodic gauge transformations 

x ,n — > 9 x c tft = c M + 9 x ,v mod 27r, c M G (-7T, 7t]. (3.6) 

Therefore, we realize the gauge fixing procedure as proposed in ||. The 
successive Lorentz gauge iteration steps are always followed by non-periodic gauge 
transformations suppressing the ZMM. At the end we check, whether the gauge 
field contains yet DDS. The latter can be excluded simply by repeating the same 
algorithm starting again with a random gauge transformation applied to the 
same gauge field configuration. We call the combined procedure zero-momentum 
Lorentz gauge (ZML gauge). It provides with 99.99% probability the global maxi- 
mum of the gauge functional. The photon propagator does perfectly agree with 
the expected perturbative result throughout the Coulomb phase. 



4 Results 

We consider qQED within the Coulomb phase at (3 values between 2 and 10 for 
k values not too close to k c . Monte Carlo simulations were carried out with a 
filter heat bath method. In order to extract the pure ZMM effect, we first apply 
the standard Lorentz gauge procedure modified by initial random gauges in order 
to suppress DDS. Let us abbreviate the notation for this modified Lorentz gauge 
procedure by LG. We compare the result with that for the ZML gauge described 
above. 

For both these gauges we have computed the averaged fermion correlator 
as defined in Eqs. ( |2.3| , |2.4|) and normalized to unity at r = 1. For invert- 



ing the Wilson-fermion matrix we employed the conjugate gradient method and 
point-like sources. In the upper part of Fig. [I] we have plotted the vector part 
< T v (t,9) > e for = 2, K = 0.122 and lattice size 12 x 6 3 . The situ- 
ation seen is typical for a wide range of parameter values within the Coulomb 
phase. Obviously, there is a strong dependence of the fermion propagator on the 
gauge copies differing by the different account of ZMM. If ZMM are present the 
propagator decays much stronger than when they become suppressed. 

The masses to be extracted seem to have different values. However, let us try 
to extract the fermion mass as it is usually done with an effective mass m e //(r) 
determined from the vector part of the free fermion propagator (Eq. (|2.7|) with 
M = ). I.e. we put 
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(IY(T + l;e)) e = cosh [ln(m e// (r) + l)(iV f /2-r-l)] 
(Mr;0)) fl cosh [\n(m eff (r) + l)(N t /2-r)] ' { ' ) 

In the lower part of Fig. [I] the corresponding numerical results for the effective 
masses are shown. In the standard LG case no real plateau is visible, whereas the 
ZML case provides a very stable one. Thus, the ZML gauge yields a reliable mass 
estimate, whereas the standard case fails. Naively, when only considering the LG 
method, one would be tempted to relate a 'bad plateau' to finite-size effects and 
to believe that the given LG effective mass result is already near to the real mass. 
Such a point of view - met in the literature - obviously fails. Taking now the 
ZML mass estimate as the reliable one the LG estimate fails by a factor ~ 3, 
in our case. 

One might ask, whether the troubles with the LG method disappear, when we 
increase (3 and/or the lattice size. In order to answer this question, we first check, 
how the ZMM-distribution changes with (3 and with the lattice size. We have 
measured the distributions of the moduli of the space- and time-like ZMM P(\(j) s \) 
and P(|</>4|), respectively, for the LG case (with DDS suppressed). In Fig. ||the 
corresponding space-like ZMM distributions are drawn. The distributions turn 
out always to be approximately step-like and bound by < (f) max ~ 7r/A^ 

with an average value 

<W> ~ ^ . (4.2) 

In order to estimate roughly the effect of the ZMM on the fermion propagator 
for varying f3 and lattice size we consider the zero-momentum approximation as 



follows. According to Eqs. (|2.6| , 2.7 ) we compute the fermion propagator only 
within the constant background modes extracted from the quantum gauge fields 
in the LG case with the distribution -P(0 M ) • Therefore, we compute 

( ' {T:n)) " <r(r;0),= /[d0]P(0)r(r;0)/ Ad0]P(0), (4.3) 



(r(l;0)V 

where T stands for the scalar or vector part of the fermion correlator ( |2.6| , |2.7|) . 



The fermion mass m is related to k according to Eq. ( |2.S| ). The results of 
this calculation for the vector part of the fermion propagator in the LG case are 
presented in Fig. |] together with the corresponding free (i.e. zero-background) 
propagator (dashed lines). One can see that the effect of the ZMM does not 
weaken with increasing (3 and lattice size, respectively. Given the average (l^l) 

as in Eq.(@3) one finds from Eqs. (£01 HI EH) that the ZMM effect does not 
disappear even in the limit — > oo . 

Preliminary Monte Carlo computations of the fermion propagator within the 
full gauge field background confirm these observations. 

We can take the approximation described above in order to check, how the 
corresponding effective fermion mass would behave. This result is shown in Fig. 
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|4j. We clearly see, that for the LG case providing the ZMM background field 
configurations we do not find a plateau (full lines). The effective mass values 
strongly differ from the real ones, i.e. mo of the free propagator (dashed lines). 

Finally, in Fig. |5] we present the fermion mass extracted from the vector 
fermion propagator within the ZML gauge for (3 = 2.0 and various K-values. We 
see a nice linear behaviour from which by extrapolating to zero mass (solid line) 
we estimate the critical value k c = 0.1307 ± 0.0001. 



5 Conclusions 

We have studied the special effect of the zero-momentum modes of the gauge 
field on the gauge dependent fermion correlator. 

We have convinced ourselves that the standard Lorentz gauge fixing prescrip- 
tion to maximize the functional (|3.2|) provides gauge copies with ZMM (besides 
DDS). These modes disturb the fermion correlator in comparison with perturba- 
tion theory and consequently spoil the (effective) mass estimate. A Lorentz gauge 
employing non-periodic gauge transformations in order to suppress the ZMM - 
additionally to DDS - (the ZML gauge) allows to reach the global maximum of 
the Lorentz gauge functional. Furthermore, it provides a reliable fermion mass 
determination, at least, if k is chosen not too close to the chiral critical line k c (/3). 
A computation of the fermion propagator with constant background gauge fields 
taken from the ZMM of the quantum fields demonstrates the disturbing effect of 
these modes very clearly. Moreover, it shows the effect to be independent of the 
bare coupling and not to disappear for large volumes. 

So far, we have studied the quenched approximation of U(l) lattice gauge 
theory. The gauge action Q2.1|) is invariant under non-periodic gauge transfor- 
mations (|3.6|). Thus, we are allowed to use the ZML gauge for evaluating gauge 
dependent objects. Contrary to the gauge action, the fermionic part (|2.2j ) does 
depend on the ZMM because of the (anti-) periodic boundary conditions. In 
this case another solution of dealing with the Gribov problem has to be searched 
for. One way, nevertheless, could be considering standard LG and taking the 
constant background modes properly into account in describing the perturba- 
tive finite-volume fermion propagator and then identifying correspondingly the 



renormalized fermion mass (see [16]). This is under consideration now. 
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Figure 1: The fermionic vector propagator (a) and the effective mass (b) at 
(3 = 2. and k = 0.122 on a 12 x 6 3 lattice for LG and ZML gauges as explained 
in the text. 
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Figure 2: Distributions of the spacelike zero-momentum mode at different (3- 
values and lattice sizes. 
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Figure 3: Free fermionic vector propagator (dashed line) and averaged constant- 
mode propagator (full line) for two /3-values and lattice sizes 12 x 6 3 , 16 x 8 3 . 
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Figure 4: Effective masses corresponding to the fermion propagator results shown 
in Fig. |3|. 
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Figure 5: Fermion mass as a function of inverse k obtained within the ZML gauge 
for (3 = 2.0 on a 12 x 6 3 lattice. The solid line represents a linear fit providing 
k c ((3) = 0.1307 ±0.0001 . 
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